b — > dds transition and constraints on new physics in B decays 
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The b — » dds transition gives extremely small branching ratios within the standard model, thus 
providing an appropriate ground for testing new physics. Using renormalization group technique we 
determine the Wilson coefficients and the mixing of the operators which contribute to the 6 — > dds 
transition. We consider contributions to this decay mode from the supersymmetric standard model 
with and without 7£-parity, as well as from a model with an additional neutral Z' gauge boson. Using 
Belle and BaBar upper bounds for the B~ — > -k~-k~ K + branching ratio we constrain contributions of 
these new physics scenarios. Then we calculate branching ratios for two- and three-body nonleptonic 
B~ meson decays driven by the 6 — > dds transition, which might be experimentally accessible. 
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I. INTRODUCTION 

Among many ongoing searches for physics beyond the 
standard model (SM), rare B meson decays seem to of- 
fer good opportunities for discovering new physics. In 
particular the experimental results on decay rates and 
the parameters describing CP violation in the B me- 
son nonleptonic two-body weak decays such as B — > ttK 
and B — > (f>Kg have attracted a lot of attention during 
the last few years (see e.g. and references therein). 
In the theoretical explanation of these decay rates and 
CP violating parameters it is usually assumed that an 
interplay of the SM contributions and new physics oc- 
curs. On the other hand, there are processes of the type 
b — > ssd and b — ► dds which are extremely rare within 
the SM. A careful study of the b — > ssd transition has 
been done 

H Hll II and the decay B~ ik + K~K- 
has been suggested as the most appropriate mode among 
possible candidates for experimental searches. The upper 
limit was first determined in Jul and subsequently con- 
strained by both B factories @,|a|. These upper bounds 
gave an unique opportunity to determine constraints on 
a variety of scenarios of new physics such as the mini- 
mal supersymmetric standard model (MSSM) with and 
without 7?.-parity violation (RPV), variations of the two 
Higgs doublets model (THDM), and models with addi- 
tional neutral gauge bosons. Using constraints from this 
decay rate the AS = 2 two body decays of B~ were 
considered 3 as well as AS — 2 decays of B c [Toj . 

The b — > dds transition has not been subject of such 
intensive theoretical studies although experimental infor- 
mation on the upper bound for the B~ — > tt~tt~K + 
decay rate already exists. Namely, the BaBar collabo- 
ration has reported that BR(B~ — > 7r~7r~A' + ) < 1.8 x 
10" 6 [3, while the Belle collaboration found BR(B~ -> 
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7r-7r-A"+) < 4.5 x 1(T 6 0. Hopefully soon the LHC-b 
would give even better constraints. 

Some time ago Grossman et al. |l2j have investigated 
the decay mechanisms of B — > A-7T decays and found 
that new physics might give important contributions to 
the relevant observables. Within their study of penguin 
operators which could receive contributions due to new 
physics, these authors also included the effects of the 
AS = — 1 transition. In their search for the explana- 
tion of the B — > Ktt puzzle, the authors of have 
investigated the B — > Ktt decay mode within a model 
with an extra flavor changing Z' boson, making predic- 
tions for the CP violating asymmetries in these decays. 
Z' mediated penguin operators have also been consid- 
ered in many other scenarios. Contributions of super- 
symmetric models with and without RPV in the same 
decay channel were discussed in Ref. 0. The difficulty 
with this decay mode is that the SM contribution is the 
dominant one. The use of quantum chromodynamics in 
the treatment of the weak hadronic B meson decays is 
not a straightforward procedure. Numerous theoretical 
studies have been attempted to obtain the most appropri- 
ate framework to describe nonleptonic B meson decays 
to two light meson states. But even the most sophisti- 
cated approaches such as QCD factorization (BBNS and 

scet) [ii m m m m nini m m minm stm 

have parameters which are difficult to obtain from "first 
principles" . Consequently, searches for new physics in 
decay modes dominated by SM contributions suffer from 
large uncertainties. 

In this paper we suggest to search for the effects of new 
physics in rare decays for which the SM gives negligible 
contributions. We only consider B~ meson decays driven 
by the b — > dds transition, since the u anti-quark is a 
spectator in this process and one should not worry about 
possible contributions of the SM penguins. The measure- 
ment of decay rates for the modes in which the "exotic" 
b — ► dds transition occurs might give an unique oppor- 
tunity to constrain parameters describing new physics. 
These constraints may then be compared with those ob- 
tained from other processes such as K° — K° and B° — B° 
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transitions. 

In the second section, we describe the b — ► dds decay 
and consider contributions of various new physics mod- 
els. First we determine the Wilson coefficients of the 
hadronic operators contributing to the effective Hamilto- 
nian in an extended operator basis, which is applicable 
for variety of the new physics scenarios. Namely, we in- 
vestigate inclusive b — > dds within the MSSM with and 
without RPV, and within an extension of the SM where 
an additional flavor changing Z' neutral boson appears. 
In Sec. Ill we write explicit expressions for the tran- 
sition matrix elements entering in exclusive nonlcptonic 
decay rates. Then we consider possible candidates for 
the experimental searches. First we study the three-body 
decay B~ — > tt~tt~ K + , which has been already investi- 
gated by both B meson factories. Then we derive decay 
rates for two-body decays B~ — > ir~K , B~ — > p~K°, 
B~ -> n-K*°, B~ -> p~K*°, and three-body decay 
B~ — * 7r~ D~ Df . In Sec. IV we comment on possibili- 
ties to observe effects of new physics in considered decays 
and summarize our results. 



II. INCLUSIVE PROCESSES 

The effective weak Hamiltonian encompassing the b — > 
dds process has been introduced by the authors of |l2j | in 
the case of B — ► Kir decays. Following their notation we 
write it as 
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where Cj and Ci denote effective Wilson coefficients mul- 
tiplying the complete operator basis of all the four-quark 
operators which can contribute to the process b — ► dds. 
We choose 
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plus additional operators ©1,2,3,4,5, with the chirality ex- 
changes L <-*• R. In these expressions, the superscripts i, j 
are SU(3) color indices. All other operators with the cor- 
rect Lorentz and color structure can be related to these 
by operator identities and Ficrz rearrangements. Wc per- 
form our calculations of inclusive and exclusive decays at 
the scale of the b quark mass (fx = mj,), therefore we have 
to take into account the renormalization group running 
of these operators from the interaction scale A. At lead- 
ing log order in the strong coupling, the operators ©1,2 
mix with the anomalous dimension matrix 
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The same holds for operators C 4 , 5 (■y(0 1 2 ) = -y(0 4 5 )) 
due to Fierz identities, while the operator O3 has anoma- 
lous dimension 7(03) = a s /ir. Anomalous matrices for 
chirally flipped operators ©1,2.3,4.5 are identical to these. 

Within the SM only the operator O3 contributes to the 
b — ► dds transition at one loop with the Wilson coefficient 
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where the functions f(x) and g{x,y) were given ex- 
plicitly in Using numerical values of the relevant 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 
from PDG |23 and including the Vtd phase one finds 
|Cf M | < 2.5 x 1CT 13 GcV" 2 . Renormalization group 
running from the weak interaction scale to the bottom 
quark mass scale, due to anomalous dimension of the op- 
erator 3 , induces only a small correction factor which 
can be safely neglected. The inclusive b — > dds decay 
width within the SM is then lil 



in',: 



\r<SM 
t sm _ l°3 

inc - 48(2tt) 3 



(5) 



which leads to the branching ratio of the order 1CU 14 . 

Next we discuss contributions of several models con- 
taining physics beyond the SM: the MSSM with and 
without RPV and a model with an extra Z' boson. For 
the THDM on the other hand, the contributions to the 
b — > dds transition coming from charged Higgs box dia- 
grams were found to be negligible. Namely, due to the 
CKM matrix elements suppression they would be even 
smaller than those found in [4| for the analogue case of 
b — > ssd. Consequently wc choose to neglect them. In 
addition, the tree level neutral Higgs exchange amplitude 
is proportional to \£,db£,ds\ / m-fi , where £db and £ds are fla- 
vor changing Yukawa couplings and m# is a common 
Higgs mass scale. This ratio is constrained from the neu- 
tral meson mixing [a. Using presently known values of 
Attik and Amg |27j one can obtain an upper bound of 
\£,db£,ds\ / m 2 H < 10~ 13 GeV -2 rendering also this contri- 
bution negligible |45j . 

In the MSSM, like in the SM, the main contribution 
comes from the O3 operator, while the corresponding 
Wilson coefficient is here 
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[24x/ 6 (x)+66/ 6 (x)], 
(6) 

as found in analyses |28( taking into account only contri- 
butions from the left-handed squarks in the loop. The re- 



cent limits on (?2 1 ^13 |29t i30l l31| disallow significant con- 
tributions from the mixed and the right-handed squark 
mass insertion terms. Therefore, we only include the 
dominant contributions given in the above expression. 
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We follow Ref. [3(| and take x = ra-gi / toj 2 = 1 and 
the corresponding values of li^is) LL ( x — 1)| < 0.14 and 
\( S 2i) LL ( x = 1)| < 0.042 [2g. We take for the aver- 
age mass of squarks toj = 500 GeV and for the strong 
coupling constant as = 0.12, and find |C| /S,S ' A/ | < 
1.6 x 10~ 12 GcV~ 2 . Using Eq. © and substituting for 
the correct Wilson coefficient one finds the MSSM pre- 
diction for the inclusive b — > dds decay branching ratio 
of the order of 10 -12 . 

If RPV interactions are included in the MSSM, the 
part of the superpotential which becomes relevant here 
is W = X' i j k LiQjdk, where and k are family indices, 
and L, Q and d are supcrficlds for the lcpton doublet, the 
quark doublet, and the down- type quark singlet, respec- 
tively. The tree level effective Hamiltonian receives con- 
tributions from the operators O4 and O4 with the Wilson 
coefficients defined at the interaction scale A ~ ms 
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The renormalization group running of the operators 
induces a common correction factor for C PPV (p) = 
fQCD(riC* pv and C* PV {p) = f QCD (p)C ppv : 
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which evaluates to /Qcn( TO fc) — 2 for a range of sneu- 
trino masses between 100 GeV < < 1 TeV. In addi- 
tion, the mixing with the operators O5 and O5 induces a 
small contribution to the Wilson coefficients C5 (fi) = 
Iqcd^)C PPV and C PPV {p) = f Q CD^)C PPV -- 
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which is of the order fQCD(mb) — 0.4 for the chosen 
sneutrino mass range. The relevant part of the effective 
Hamiltonian we use in this scenario is then 
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theories with large extra dimensions |33l |. This induces 
contributions to the effective tree level Hamiltonian from 
the operators as well as ©1,3. Following H2,|33, the 
Wilson coefficients for the corresponding operators read 
at the interaction scale A ~ mz> 



We neglect the /qcd suppressed contributions of O5, O5 
to the amplitudes in the cases where the operators O4, 
O4 give non-zero contribution. The inclusive b — > dds 
decay rate induced by the RPV model becomes 
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Present experimental bounds on the individual RPV cou- 
plings contributing to the effective Wilson coefficients 
C PPV and C PPV do not constrain this mode, and we 
extract the bounds on the relevant combination from ex- 
clusive decays in Sec. IIVI 

In many extensions of the SM [33] an additional neu- 
tral gauge boson appears. Heavy neutral bosons are also 
present in grand unified theories, supcrstring theories and 
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where y — (g2/gi) 2 (pi sin 2 9 + P2 cos 2 9) and pi — 
myp/ m i cos2 @w ■ In this expression gi, g2, mi and m,2 
stand for the gauge couplings and masses of the Z and Z' 
bosons, respectively, while 9 is their mixing angle. Again 
renormalization group running induces corrections and 
mixing between the operators. As already mentioned, the 
mixing of operators ©1,2 and their chirally flipped coun- 
terparts is identical to that of operators ©4,5 since these 
operators arc connected via Fierz rearrangement. Thus 
the same scaling and mixing factors /qcd and Jqcd ap- 
ply. For the operator O3 on the other hand the renor- 

(M) = /qcd(m)C 3 z ' with 



malization can be written as 
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(13) 

In particular for a common Z' boson scale of mz> — 
500 GeV one gets numerically /qcd( to 6) — 2, 
fQCD{m b ) ~ 0.4 and /g CD (m 6 ) ~ 0.8. The full con- 
tributing part of the effective Hamiltonian in this case 
is 



n_j Z 



/qcd(m) [cf'Oi(/i) + <5f'Oi(/x)" 
+/qci?(m) [Cf ' 2 (m) + Cf 2 (m) 
+/q CO (m) Icfo^ + Ci'd^)} . (14) 



For the inclusive b — > dds decay rate the 02 and 02 are 
numerically suppressed due to the Jqcd factor and we 
write 
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In Sec. HVI wc discuss bounds on Wilson coefficients Cf 3 



and C*f 3 which might be estimated from the B 
Tt~7r~ K + decay rate. 



III. EXCLUSIVE B~ DECAY MODES 



In calculating decay rates of various B meson decay 
modes based on the b — > dds quark transition, one has to 
calculate matrix elements of the effective Hamiltonian op- 
erators between meson states. As a first approximation, 
we use the naive factorization of three-body amplitudes, 
and express the resulting two-body transition amplitudes 
between mesons in terms of the standard weak transi- 
tion form factors (|A1IA3J) . as dictated by the Lorentz 
covariancc. For the B — > 7r(p) transitions we use form 
factors calculated in the relativistic constituent quark 
model, with numerical input from the lattice QCD at 
high momentum transfer squared For the D s — > D 
and K — ► it transition form factors we use results of 
Refs. 0, |3j| where heavy meson effective theory and 
chiral Lagrangian approach were used. 

For the decays of the B~ meson to three pscudoscalar 
mesons P, P\ and Pi we first derive a general expression 
for the factorized matrix element of the O3 operator, rel- 
evant in the framework of SM (MSSM) 



(P 2 (K)A(pi)|<M°) (P(p)\dY l b\B-(p B )) = (t- u)F p ^(s)F PB (s) 

*k){rn%-m a P ) [ f p 2 p 1{s)f pb {s) 



(m 2 Pi 



F P ^(s)F PB (s) 



(16) 



Because only vector currents contribute in the above 
expression, it also applies for operators O3, 0\ and 
0\. Form factors F\ and Fq are defined in the Ap- 
pendix A and the Mandclstam kincmatical variables are 
s = (pb - p) 2 , t = (p B - Pi) 2 and u = (p B - p 2 ) 2 . 

In the context of RPV, the contributions of the oper- 
ators 04 and 04 to hadronic amplitudes are dominant. 
One can use the Dirac equation to express scalar (pseu- 
doscalar) density operators in terms of derivatives of vec- 
tor (axial-vector) currents 
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(17b) 



Using these relations we derive an expression for the fac- 
torized matrix element of the 04 and 04 operators, con- 



tributing only with their scalar parts 
(P 2 (P2)Pi{pi)\ds\0)(P(p)\db\B-(p B )) 



(m 2 Pi 



'P-2 



){m 



-^F p ^(s)F™(s). (18) 
(m b - md)(m s - m d ) 

In the case of the Z' model one encounters contribu- 
tions of the operators 0i,2,3 and 0i,2,3- The color non- 
singlet operators 02 and 2 can be Fierz rearranged to 
04 and 04 and then Eq. (|18fl applies as well. Remaining 
operators are all of the V±A form and their contribution 
to the amplitude is already given by Eq. (|16f) . 

In two-body decays with a vector meson V and a pseu- 
doscalar meson P in the final state we sum over the po- 
larizations of V. The sum in our case reduces to 
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where ey is the polarization vector of V and A is defined 
as \(x, y, z) = (x + y + z) 2 — A(xy + yz + zx). 
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For decay to two vector mesons in the final state we 
use the helicity amplitudes formalism as described in 
Ref. Non-polarized decay rate is expressed as an 

incoherent sum of helicity amplitudes 



r = -pi T (\H a \ 2 + \H +1 \ 2 + \H_ 1 \ 2 

8irm B \ 



(20) 



where pi is momentum of the vector meson in B meson 
rest frame and helicity amplitudes arc expressed as 



y/\(m 2 g ,ml,m$) 
±i±i = a ± c, 
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Vector meson masses are denoted by mi^, while defini- 
tion of the constants a, b and c is given by general Lorentz 
decomposition of the polarized amplitude 



Hx = eUA) e L(A) a«T + 
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where ei^ and 77^2 are the vector mesons polarizations 
and momenta. 



Hadronic matrix element entering in the amplitude 
for B~ -> tt~tt~K + in SM (MSSM) is readily given by 
Eq. (JTBJ) after identifying P = tt~ , P 1 = K+ , P 2 = 7T - 
and using appropriate form factors given in the Ap- 
pendix A. Eq. H18f) is used instead for RPV, while the Z' 
amplitude incorporates both Eqs. lfR)|> and JTSJ . There 
are two contributions in each model to this mode, with an 
additional term with the u <-> s replacement in Eqs. i|16|) 
and 118|) , representing an interchange of the two pions in 
the final state. After phase space integration, the decay 
rates can be written very compactly with only Wilson 
coefficients left in symbolic form: 



B. B 



7v-D~Dt 



In calculation of the B~ — ► tt~ D~ Df decay rate again 
we use Eqs. (I16[) and l|18|) now with substitutions P = ir~ , 
Pi = and P2 = D~ . Numerically this yields 
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1-5 



(26) 
(27) 



x 1.5 x KT° GeV° 
2 x 5.5 x 1CT 9 GeV 5 



Re 



(cf+Cf) (cf +G 



x 5.7 x 10~ 7 GeV 5 . 



(28) 



These decay rates are suppressed due to the small phase 
space in comparison to the rates of the B~ — > ir~ir~ K + 
decay. 



7T~K 







The authors of Ref. 0] addressed this decay mode as 
the wrong kaon mode, being highly suppressed in the SM 
compared to the decay with K° in the final state. The 
operators (D\^ and Oi^ that are present in SM (MSSM) 
and Z' model have the following contribution: 

(K°(p K )\d 7fi j 5 s\0) (v-(p*)\drfb\B-(j>B)) 
= i(m B -ml)f K FZ B (m 2 K ). (29) 

Operators O4 and O4, relevant for the RPV and Z' mod- 
els result in 

(K°(p K )\d 1 5 S \0)(7r-(p^\db\B-(p B )) 
im 2 K (m% - ml) 



(m b - m d )(m s + m d ) 



(30) 



However, in the latter two models, the two chirally 
flipped contributions to the amplitude have opposite 
signs, resulting in a slightly different combination of Wil- 
son coefficients in comparison with the B~ — > 7r~7r~ K + 
decay rate 
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D. B~ 



P-K° 



Using form factors parameterization (|A3|) of the pseu- 
doscalar to vector meson transition we derive the follow- 
ing two factorized matrix elements of axial-vector and 
pscudoscalar operators: 

(K°(p K )\d 7 ^ 5 s\0) (p-(e p ,p p )\dn 5 b\B-(p B )) 

= -2m p f K Af(m 2 K )e* p - PB , (34) 

(X (p A -)|d7 5 .s|0) <p-(e p ,p p )|d 7 5 6|i3-(p B )> 

2m p m 2 K 



-t/kAP (m K )e*-p B . (35) 
{mb + md)(m s + md) 

Finally, we sum over polarizations of the p meson using 
Eq. I|19fl . and the unpolarized decay rates read 
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x 1.5 x 1(T 3 GeV 5 . 
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E. B~ 



TT-K*° 



Factorized matrix element is here a product of vec- 
tor meson K*° creation amplitude ljA2b(l and B~ — ► tt~ 
transition amplitude. Operators which involve vector 
currents result in 

(K*°(e K , PK )\d^ S \0) (ir-(p n )\d^b\B-(p B )) = 

2g K *F? B (m 2 K ,y K - PB , (39) 

while the density operators O4 and O4 do not contribute, 
as a result of Eqs. l(T7J) and (|A2 b|) . Consequently, in the 
RPV model this mode is dominated by the operators O5 
and 05 which are, as mentioned in Sec. ^ suppressed 
by the renormalization group running. Using Ficrz rear- 
rangements, we write them down as <D\, 0\ and yield an 
additional 1/2 suppression factor. 
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x 4.7 x KT 4 GeV 5 
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2.3 x 1CT 3 GeV 5 . 



(42) 



Like in the previous case, this mode only receives con- 
tributions from the renormalization group suppressed 
RPV terms. Wc calculate unpolarized hadronic ampli- 
tudes of the operators Oi,3 and 01,3 by utilizing the he- 
licity amplitudes formalism. Using form factor decompo- 
sition (|A2bl IA3(I , we write down the expression for the 
polarized amplitude (|22|l and identify constants a, b and 



b= - 



(m B + m p )g K , Af(m 2 K .)(C-C), (43a) 
g K ,Af(m 2 K *)(C-C), (43b) 



4 v-~ ' ' 
1 m,K*m p 



2 m B + m f 

c = -\-^^m^ pB {rn 2 K .){C + C). (43c) 
l m B + m p 

C and C are combinations of the Wilson coefficients 

f AI S) S M 

present in a considered model. Wc have C = C3 , 
C = in the SM (MSSM), C = -f Q CD(m b )C£ PV /2, 
C = -f QC D(m b )C^ PV /2 in the case of the RPV 
model and C = fQCD{m b )Cf + ]' QCD (m b )C% ' , C = 

fQCD(m b )Cf + f' QCD {m b )Cf in the Z' model. Decay 
rates are then 



^(MS)SM 
pK* 



c 



(MS)SM 



^RPV 
pK* 



T Z' . 
1 pK* — 



c PPV 



qRPV 



x 9.2 x 10" 

2 



GeV 5 



(44) 



RPV 



C 



cf + cf 



c 



RPV 



x 1.4 x 1(T 6 GeV 5 
2 

1-5 



x 3.5 x HP 5 GeV , (45) 



Cf -Cf 



Ci 



Co 



x 1.4 x KT 4 GeV 5 



x 3.5 x KT 3 GeV 5 



x 2.2 x 1CT 5 GeV 



cf - cf 



x 5.7 x 1CT 4 GeV 



Re 



Cf + Cf 



x 1.1 x KT 4 GcV d 



Re 



Cf -Cf 



Ci 



cf 



Ci 



x 2.8 x KT 3 GeV 



IV. DISCUSSION AND RESULTS 



(46) 



We have investigated the b — > dds transition within 
the SM, MSSM without and with RPV terms and within 
a model with an extra Z' gauge boson. The SM contri- 
bution leads to extremely small branching ratio for this 
transition. 

First we have calculated the effects of QCD on the Wil- 
son coefficients caused by the renormalization group run- 
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ning. The moderate increase of the MSSM compared to 
the SM predictions is still too insignificant for any exper- 
imental search. The MSSM with RPV terms, however, 
might give significant contributions and a possibility to 
shrink down the parameter space even further. The Z' 
model exhibits its structure through interplay of different 
interaction scale couplings and might also give opportu- 
nity to constrain its relevant parameters. In the case of 
the two Higgs doublet model we do not expect any sizable 
effect as already noticed in the case of b — ► ss d decays Q. 

In the b — > dds decay a particular combination of the 
model parameters appear which can be constrained us- 
ing the B~ — > 7r~7r~A' + decay mode. In our calculation 
we have relied on the naive factorization approximation, 
which is as a first approximation sufficient to obtain cor- 
rect gross features of new physics effects. One might 
think that the nonfactorizablc contributions might in- 
duce large additional uncertainties, but we do not expect 
them to change the order of magnitude of our predic- 
tions. Additional uncertainties might originate in the 
poor knowledge of the input parameters such as form 
factors. However, we do not expect these to exceed more 
than 30%. 

Using the stringest experimental bound for the 
BR(£?~ — ► -k~tt~K + ) < 1.8 x 10~ 6 and normalizing the 
masses of sneutrinos to a common mass scale of 100 GeV 
we derive bounds on the RPV terms given in Eq. J7J 



'100 GeV N 2 



E 



( A 'n31 A 'nl2 + A 'n21 A nl3) 



< 8.9 X TO" 5 . 



(47) 



Complementary bounds coming from measurements of 
K° — K° and B° — B° mixings have been established in 
Refs. 01 



Re 



Im 



Re 



^ 100 GeV \ 2 

2 „ 1 ~ J A n31 A nl2 



.n—1 v 

3 

E 



'100 GeV x 



.71=1 

3 

E 



/100 GeV 



A n31 A nl2 



A 'n21 A 



nl3 



< 2.6 x 10" 6 , (48a) 

< 2.9 x 10~ 8 , (48b) 

< 2.9 x 10~ 4 . (48c) 



From Eqs. (|48a|l and l|48b|) it becomes apparent that the 
A r n31 \'* 12 term is negligible in Eq. (|4*7|) . and the bound 
becomes simpler 



A f WO GeV \ 2 

/ , A n21 A nl3 



< 8.9 x 10 



-5 



(49) 



now being more restrictive than Eq. Q48cj l. obtained from 
B° - B° mixing. 

Assuming that new physics arises due to an extra Z' 
gauge boson we derive bounds on the parameters given 
in Eq. i|12|) • Making the simplest assumption, we neglect 
interference between Wilson coefficients (third term) in 



Eq. 1)25^. Experimental bound of this simplified ex 
pression now confines f|Cf + C{ 



Cf\) to lie 



within an ellipse with scmiminor and semimajor axes as 
upper limits 



B 12 B 13 



y B 



12 -°13 



B 12 B 13 



D \2 D \3 



< 2.6 x 10" 



< 7.1 x 10" 



(50a) 
(50b) 



Complementary bounds, involving the same couplings 
and y, originate from meson mass splittings and CP vio- 
lation in kaon system and have been derived in Ref. |32l | 
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V tte[(tf 12 -j-j < 10" 
y\Rc[(B d ^ L ) 2 } <6xl0" 

». T.„ U r> d R,L\2i , o m-11 



(51a) 
(51b) 
(51c) 



To combine those bounds with Eqs. I|50(l . one should ab- 
sorb dimcnsionless y into the coupling constants by redef- 
inition B^;^ = ^12(13) ■ However, when we include 
the constraints given by Eqs. (|50|l . we obtain no further 
improvement of the bounds on individual B couplings. 

Nevertheless, the bounds (|47|l and l|50|) are interesting 
since they offer an independent way of constraining the 
particular combination of the parameters, which arc not 
constrained by the B° d -B° dl B°-B°, K°-K° oscillations, 
or by B~ — > K~K~tt + decay rate (see e.g. 0). 

Using these inputs we predict the branching ratios 
for the various possible two-body decay modes and the 
B~ — > n~D~Df decay. Applying bound (|4*7j) to the 
RPV model decay rates is straightforward, except for the 
B~ — > ir~K° and B~ — > p~K*° decay modes. In order 
to make predictions for these two modes, we assume as 
in H||, that interference term C? PV C? PV * is negligi- 
ble, which leads to the approximation \C PPV — C PPV | ~ 
\C? PV +C? PV \. 

In the case of the Z' model, there are contributions 
from Wilson coefficients "1" (Cf' , Cf' ) and "3" (Cf' , 
C% ). We have already neglected the interference terms 
between "1" and "3" in Eq. (23 to obtain bounds lfSTj|) 
and we assume that these terms are small for all con- 
sidered decay modes. Using Eqs. l|50f) we can now pre- 
dict branching ratios for decay modes B~ — » tt~ D~ Df , 
B~ — > p~K°, and B~ — > ir~K . The remaining two 
decay rates B~ — > tt~K° and B~~ — > p~K*° can be con- 
sidered after we neglect interference terms Cf Cf * and 
Cf Cf * . The results are summarized in Table [I] 

The SM gives negligible contributions. The MSSM 
is increasing them by two orders of magnitude, which 
is still insufficient for the current and foreseen exper- 
imental searches. Using constraints for the particu- 
lar combination of the RPV parameters present in the 
B~ — > ir~ir~~K + decay we obtain the largest possible 
branching ratios for the two-body decays of B~ — > p~K° 



8 



Decay 






SM 




MSSM 




RPV 




Z' 




BT -> 


Tf" 


TT-K+ 


3 x 10" 


16 


1 x 10" 


14 


— 




— 




B- -> 


7T" 


D~D+ 


1 x 10" 


21 


6 x 10" 


20 


2 x 10" 


-8 


3 x 10" 


-9 


B~ -> 


Tf" 


K° 


6 x 10" 


17 


3 x 10" 


15 


1 x 10" 


7 


5 x 10" 


7 


B~ -> 


P~ 


K° 


6 x 10" 


17 


3 x 10" 


15 


1 x 10" 


7 


8 x 10" 


7 


B~ -> 


7T" 


K*° 


1 x 10" 


-16 


5 x 10" 


15 


6 x 10" 


9 


1 x 10" 


-6 


B~ -> 


P~ 


K*° 


1 x 10" 


-16 


6 x 10" 


15 


7 x 10" 


-9 


1 x 10" 


-6 



TABLE I: The branching ratios for the AS = — 1 decays of 
the B~ meson calculated within SM, MSSM and RPV models. 
The experimental upper bound for the BR(B~ — > tt~ ir~ K + ) 
< 1.8 x 10~ 6 has been used as an input parameter to fix the 
unknown combinations of the RPV terms (IV column) and 
the model with an additional Z' boson (V column). 



and B~ -> ir-R , while for the B~ -> tt^A'* and 
B~ — ► p~K*° the RPV contribution is suppressed by 
the renormalization group running. Their order of mag- 
nitude is 10~ 9 and thus still experimentally unreachable. 
However, these two decay channels are most likely to be 
observed in the model with an additional Z' boson, if we 
assume that interference terms are negligible. 

Since in the experimental measurements only Kg or 
Kl are detected and not K° or K°, it might be difficult 
to observe new physics in the B~~ — > it~ K° decay mode. 
Namely, the branching ratio BR(B~ — > ir~Kg) = (12. 1± 
0.7) x 10~ 6 is two orders of magnitude higher than 
our upper bound for the BR(B~ — > ir~ K°) making the 
extraction of new physics from this decay mode almost 
impossible. Therefore, the two-body decay modes with 
K*° in the final state seem to be better candidates for 
the experimental searches of new physics in the b — > dds 
transitions. 
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APPENDIX A: FORM FACTORS 

We use the standard form factor parameterization of 
hadronic current matrix elements |40t l4ll | between two 
pseudoscalar mesons 



(ft(P2)|arAi|^(pi)) 

= F 1 ( q 2 ) ((px+ P2 r- 

+F a (q 



2 ,m 2 Pi -m% 

9 ) 



(Al) 

where = (pi — p 2 Y- We also use the standard decay 
constants of the pseudoscalar K and vector K*° mesons 



(K°(p)\d^ 5 s\0) = i/kp", 
(K* (e K ,p)\dn 5 s\Q)=g K ,e 



K ■ 



(A2a) 
(A2b) 



Matrix element between a pseudoscalar and a vector me- 
son is decomposed as customary 



(V{e Vl p 2 )\q^qAP{ Pl )) 



mp + my 



(A3a) 



2mv 



(V(e V ,p 2 )\q J ^ 5 q l \P( Pl )) = *e y ■ q^fq»A (q 2 ) 



+ i(mp + my) 

~ l i \ 

(mp + mv) 



A!(q 2 



{P1+P2T 



m% - m 2 v 

o q 



A 2 (q 2 ). 

(A3b) 



For the B~ — > tt~ and B~ — > p~ transitions we 
use form factors calculated in the relativistic constituent 
quark model with numerical input from lattice QCD at 
high q 2 



J 



F? B (q 2 ) 
F^(q 2 ) 



F? B (0) 



(1 - q 2 /m 2 B ,)[l - aiq 2 /m 2 B 



. Fj^O) =0.29, (Ti=0.48, 



1 — <J\q 2 /m B , + a 2 q 4 jm 4 



, ^(0) = 0.29, <j\ = 0.76, cr 2 = 0.28, 



(A4a) 
(A4b) 
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V pB (Q) 

^ ) = __<^__ 5jj Af(0)=0.30, .1 = 0.54, (A5b) 

? fl ^ = i- glgV ^w/H. ' < (0) = - 26 = CT - - 73 ' - = - 10 ' (A5c) 



<V) = - 2 2 ^ 4 , Af(0)=0.24, ^ = 1.40, ^ = 0.50. (A5d) 



The transition form factors between heavy mesons D s — ► Z) have been calculated in the chiral Lagrangian approach 
by the authors in Ref. 

F 1 DD °(q 2 )=0 1 (A6a) 



.,2 



(Stt / 4 ) /if (1430) ^fmoJ^D 



m D s m D <T - (1430) + *V 9 r^(1430) 

The same method has been used to obtain the light to light K~ — ► 7r~ meson transition form factors in Ref. 

F\ K {q 2 ) = 2 9VKW2)9K' ^ (A7a) 

q 2 -m 2 K{S92) + iyj q 2 T A '(892)(<? 2 ) 

V * ' 9 9 , ■ / / 9 \ ™. 2 -™ 2 O 9 • I 9"t-i / 9 \ ^ ' 



_ 2 9VK(&92)gK* (1 ~ q / m K(892)) g 2 Jg(1430)gS-R-(1430) 

I 2 ~ m if(892) + ^ \/<F ^ if(892)(<2 ,2 ) m if -~ m lq 2 ~ ^(1430) + * V^if (1430) (? 2 ) 



Here the decay widths of the resonances A* (892) and A(1430) are taken to be energy dependent [35| 
Tk gg2 ,2s = ( m K(892) \ 1 ( [q 2 - (mK + m^) 2 } [q 2 - (m K - m^) 2 } \ ' ^ ^ 

K{892) \ ^ J \ [ TO /^(892) _ i 171 * + m ^ 2 ]\ m K(S92) ~ ( m K ~ m ^ J A (892) ' 

r , 2 , _ ( 4(1430) \ ^ ( [q 2 - (m K + m^) 2 ][g 2 - (m K - m^) 2 ] \ 1/2 

r " (1430)( ' ?) "l v ^^J 1 14(1430) - (»* + ^) 2 ][4(1430) " (m K - m.) 2 ] ) TK ^- (A8b) 



APPENDIX B: NUMERICAL PARAMETERS 



TUB 


TUB* 


m,r 


m p 


m K+ 


m K o 


m if»°(892) 


m D 


m Ds 


«1lT(1430) 


m b 


m a 


m d 


5.27 


5.32 


0.140 


0.77 


0.494 


0.498 


0.892 


1.87 


1.96 


1.41 


4.2 


0.10 


0.006 



TABLE II: Meson and quark masses (in GeV) used in our calculations are taken from PDG (271 . 



Decay constants of the pseudoscalar K° and vector K*° mesons are / A = 0.160 GeV [34[ and $k * = 0.196 GeV 2 
respectively. Further numerical parameters relevant for the D s — > D and A' — > 7r transitions are [Tol . l27j | = 3.73. 
/if (1430) = 0.05 GeV, r K(1430) = 0.29 GeV, g SK (u30) = 3.7 GeV, g VK (s92) = 4.59, and T K(sg2) = 0.051 GeV. 
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